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Abstract 

A general case of a spatially nonuniform planar layered Ising model, or an 

equivalent quantum Ising chain, is analysed with an exact functional real 

space renormalization group. Various surface, finite size, quasiperiodic and 

random layer (McCoy- Wu) universality classes are obtained and discussed 

within a single theoretical framework leading to new insights into the nature 

of random criticality. 
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While equilibrium critical phenomena in spatially uniform systems are gradually becom- 
ing a standard textbook chapter, much less is known about the way the critical state of 
matter is altered by the presence of spatial inhomogeneities of a general form. Models with 
uniformly distributed point disorder, on one hand, and with a planar surface, on the other, 
have attracted considerable attention and considerable progress in understanding of both 
has been achieved. However, the developments along these two directions have been largely 
independent, leaving a whole range of inhomogeneities largely unexplored. Important ex- 
amples falling into this range are, to name just two, disordered quantum systems in the 
path integral representation of which the disorder is fully correlated in the time direction, 
and fluids in porous media , in which case disorder enters in the form of random surfaces. 

Here I am attempting to bridge this gap by presenting a new analysis of the class of 
layered planar Ising (LPI) models. A layered Ising model is defined as having its bond 
strengths depending in an otherwise arbitrary way on just one of the d spatial coordinates, 
z; in each layer perpendicular to that coordinate the bonds are of the same strength. These 
models may be viewed as a natural step in generalization of models with one planar surface 
towards the complexity of real systems. In particular, any LPI model is equivalent to an 
s = 1/2-quantum Ising chain in transverse field which, depending on the way the transfer 
matrix is applied, can be made inhomogeneous either in space or in time. Besides, the 
rapidly developing technology of artificially grown multilayers provides natural experimental 
realizations of the layered geometry. Starting from the pioneering work of McCoy and Wu 
a number of LPI models, both with surfaces and randomness 0, have been solved 
exactly. Subsequent scaling analysis of the exact solutions for surface and film geometries 
has played a crucial role in formulating surface and finite size scaling 0. More recently 
important progress has been made towards understanding of the random layer model 
so the whole class seems to be a good testing ground for unifying ideas. 

This work builds upon the microcanonical density functional method, proposed recently 
by M. E. Fisher and the author. Exact Euler-Lagrange equations obtained for the LPI class 
in P] are used to construct a functional renormalization group (RG) flow in the space of 
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energy density profiles. Solving for the fixed points of this flow generates a list of universality 
classes found in LPI models. Common features of the flow in the vicinity of the fixed points 
describing the free film and random criticalities leads to interpretation of the latter as a 
hierarchy of thin-thick film crossovers with a significant element of quasi-one-dimensional 
behavior. This qualitative insight allows then for quantitative predictions to be made re- 
garding critical behavior in quasiperiodic and long-range-correlated random sequences. 
The Euler-Lagrange equations found in 

describe the response of the energy density components eq{z) to the scaled temperature field 
t{z), which characterizes deviation from criticality of the local bond strengths in a given 
LPI model. Here q is the wave number in the direction parallel to the layers, Bq{z) is the 
contribution to the total energy density from the excitations belonging to the g-sector of the 
equivalent fermion problem |^,^; iq{z) = dsg/dz, etc. An inverse latttice constant has been 
absorbed into the definition of t, giving it the proper scaling dimension of inverse length; the 
normalization of makes the absolute value |t| = ^^^{z), where C,{z) is the correlation 



length of a uniform model with the same bond strengths as the given z-layer (recall that 
^ oc \t\^^ in a uniform planar Ising model). In this formalism the vertical bonds (connecting 
subsequent layers) have to be ferromagnetic, but the horizontal bonds can take either sign 
thus allowing for restricted forms of frustration (cf. 0). Eqs. 1 are exact in the scaling limit 
t{z) <^ A, with the cutoff A being of the order of the inverse lattice spacing. They have to 
be solved for every q between and A; after that the energy density, free energy, etc., are 
obtained by integration over q of simple functions of eq{z) and eq{z). This program has been 
explicitly carried out for simple surface and superlattice geometries, as represented 
by step-function and periodic t{z), respectively. 

In the latter case, apart from the information obtained from the explicit expressions for 
the free energy and total energy density profiles, it turned out to be instructive to study 
evolution of the partial energy density profiles eq{z) in the long wave length limit g — 0; 



that evolution resembled many fetures of a functional renormalization group. Indeed, taking 
a partial derivative of the left hand side of (1) with respect to the logarithmic length scale 
along the layers, / = ln(A/g), one obtains a partial differential equation depending only on 
£5, its first two spatial derivatives, and, generally, on q — Ae~' as well. In this approach, the 
field t can be viewed as fixing the initial data at large q — A., where s{z; I — 0) ^ t{z)/A. The 
evolution of profiles as the logarithmic length scale / increases follows then its own dynamics, 
independent of the source t, which, according to (1), becomes an integral of motion for this 
flow. The simplest illustration is provided by the uniform limit t — const, de / dz — 0. 
Differentiating both sides of (1) with respect to / gives then 



As anticipated, Eq. 2 has three fixed points: £ = 0, ±1. In the scaling regime, t -C A, the 
flow always starts in the vicinity of the ultraviolet-stable critical fixed point £ = and then 
takes the system either to the high-temperature, £ = 1, or to the low-tmeperature, e — —1, 
fixed point, depending on whether the initial value e{l = 0) = t/A is positive or negative, 
correspondingly. The eigenvalue = 1 governing the flow out of the critical flxed point, 
de/dl XfE, gives the correct value of the correlation length exponent v — l/\t — l. 

The general case of position dependent t requires supplementing flow in q with rescaling 
£ — > £ exp(— a;^/), z ^ z exp(— a;^/), where the scaling dimensions a;^, ujz are tuned to achieve 
nontrivial fixed point structure and make the RG equations /-independent. In pracitice, it 
is easier to work with the integral of motion (1): each RG fixed point profile 



has to satisfy (1) with the same t{z) for every small q. A straightforward analysis yields 
several fixed points. 

First, at a;^ = 0, = 1 one finds the profile 



de/dl = e{l-e^). 



(2) 



e,{z) = q^^S{q^^z) 



(3) 




(4) 
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describing relaxation towards the uniform critical fixed point profile e = 0. The thermal 
potential 

g{z) = 2 rt{z^)dz,, (5) 



will play a crucial role below. The condition of criticality in a layered Ising model discussed 
in can be made more precise now: for (8) to make sense there has to exist a choise of 
integration constant in (5) such that the integral in (4) is bounded for all y (note that e 
cannot exceed unity). This rather stringent constraint is satisfied in the critical state of a 



periodic |]T0[ and, as one will see below, simple quasiperiodic LPI models, but not at the 
ferromagnetic transition in the random layer model, where g{z) performs a random walk. 

If positive values of t dominate, so that aX z ±oo the potential g diverges as ±1^1, 
respectively, then the profile relaxes towards the high-temperature fixed point e = 1. Similar 
relaxation to e = —1 occurs for mostly negative t; the two can be combined into the scaling 
form e = ±1 =1= q'^£±{z) + o{q^), where 

8± = - dzi dz2 e^[^(^2)-^(^^)l. (6) 

Z J —QO J Z 

The scaling dimensions describing this leading irrelevant behavior at the noncritical fixed 
points are evidently = 0, ujg = 2. 

If t{z) is mostly positive at z > and mostly negative at z < 0, so that g has a global 
minimum around z = 0, one finds a nontrivial fixed point 

S{y) = -1 + 2 / e-^^'Mz/ / e-<^^'^dz (7) 



oo J — oo 



scaling with uo^ = oJi. = 0. It describes an interface between the low-temperature phase on 
the left and the high-temperature phase on the right; note that inversion of this solution, 
z —z^ implies inverting sign of g. In fact, by choosing finite instead of infinite limits of 
the integrals in (6), (7) one still obtains a local solution to (1); those local solutions will be 
used below in analyzing the case of a generic, sign indefinite t{z). 

The remaining fixed point profiles describe response of a critical system to sur- 
face/interface perturbations. The choise of tu^ = 1/a > 1, = yields a solution 



8{y) = {I + A-YT''' (8) 

describing local, eg{z) = t{z)/{t^{z) + g^)^/^, response to a power-law tail t{z) = Az~°'. 
Scaling in this case is somewhat spurious: the local dependence of Eg on t is realized whenever 
the source t varies slowly on the scale of the local correlation length ^{z) = t^^{z), i. e. 
t^'^t <^ 1. This condition, bound to fail near every zero of t, yields a < 1 for power-law 
sources t oc z~"'. Physically it implies that the essential correlations decay faster than the 
bond strength changes, therefore any physical property can be calculated by simply adding 
contributions from quasihomogeneous regions. The marginal value a = 1 leads one to the 
scaling dimensions cj^ = 1, cj^ = of the uniform critical model. Although an explicit 
solution of the fixed-point equation 

[S-\£][l-S' + lsr'^'=q-H{y/q) (9) 

for q~^t{y/q) = A/yis not available at the moment, presumably such solution exists for any 
value of the dimensionless amplitude A\ interesting properties of this line of fixed points are 
discussed in [0. Besides t oc 2;"^, the right hand side of (9) is g-invariant for t{z) = and 
t{z) = goS{z); the corresponding fixed point profiles are 

£s = const e^'^y; £d = tanh(^o)e"^'^'. (10) 

The first one describes a surface of a critical half-plane in the absence of long range per- 
turbation in the bond strengths, including t oc z'"" with a > 1 [0. Adding two of these 
half-plane solutions yields £d describing response of a critical system to a defect line 0. 

These fixed points allow for a compact descritption of a general case of an interface 
geometry 0. In the most generic case the interface separates two noncritical half-spaces 
characterized by temperature fields ^1,^2- Only the relative sign of the ti and t2 turns out 
to matter then: ti and t2 of the same sign inevitably lead to profiles flowing into one of the 
uniform fixed points (6); opposite signs of the temperature fields, on the other hand, lead 
to the flow into the nontrivial fixed point (7). This discontinuous change in the form of the 
profile as one of the half-planes passes its critical point is at the origin of the surface latent 
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specific heat found in the exact solution of the half- infinite model @^J^. At criticality of 
one of the half-planes the flow can result in one of three fixed points, determined by the 
range of the surface perturbation in the t-field, as explained above and in agreement with the 
analysis in [j^. Finally, when both half-planes go critical simultaneously the total strength of 
the bond perturbation in the interfacial reion, as encapsulated in the amplitude Qq = 2 J tdz 
of the (5-function, determines which of the line of fixed points Sd (10) attracts the flow. 

An insight into the nature of some of these fixed-point profiles is provided by the exact 
equivalence between each g-sector of an LPI model and a ficticious one-dimensional 
classical Ising chain at temperature Td=i ~ | lng|^^ and nonuniform reduced magnetic field 
Hd=i{z)/kBTd=i ~ The energy density components Eg map onto the magnetization 

density of the corresponding g-chains. The low-g limit of the planar model corresponds 
to the low-temperature limit of the chain. The latter is known to be determined by the 
statistics of a dilute gas of pairs of antiparallel spins which I will call kinks below. A "spin 
down-spin up" kink, to which orientation a = 1 will be assigned, has to be followed by an 
"up-down" one, a = —1, etc. The energy of a kink at position z along the g-chain is a 
sum of the cost of creating the kink, ln(l/g), and of interaction with magnetic field, CFg{z), 
leading to the Boltzman factor gexp(— crfyf). Here the additive constant in the definition 
of the potential g is arbitrary but has to be the same for all kinks. One can see now 
that (7) represents a single kink localized near a global extremum of g. The uniform fixed 
points naturally correspond to the chain either fully magnetized in one of the two possible 
directions, e = ±1, or maintaining zero average magnetization, e = 0, at Td=i = 0. The 
low-g corrections to these states represent a background of thermally excited single kinks in 
(4), and pairs of kinks in (6). Other fixed points do not have so simple interpretation, as 
kinks strongly interact with each other there. 

None of the fixed points considered so far addresses directly the situation in which t{z) 
has more than one zero, so that g{z) has many maxima and minima. A crucial insight 
is provided by consideration of the simplest system with two boundaries. In the present 
formalism a strip of width L is represented by t{z) = ti at \z\ < L/2, surrounded by 
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t = +00 background [3,|T3]. At t > the energy profile is asymptotically fiat, but at ti < 
and for sufiicient thick films, \ti\L ^ 1, the contradiction between the tendency to ordering 
within the film and the absence of any order outside of it is resolved by creation of a pair of 
kinks separating "cool" interior from the "hot" background. Formally a solution is obtained 
by matching two kink solutions (7) of opposite orientations with three noncritical solutions 
(6) by choosing appropriate integration constants in those. This construction collapses, 
however, when q is decreased past ^|~^ oc e'*^'^. In the kink picture this comes naturally: the 
kink-antikink pair is annihilated as soon as the "self-energy" of the pair 2| lng| exceeds the 
potential difference \g{L/2) — g{—L/2)\ = 2\ti\L supporting the pair. Physically ||10|, the 



exponential dependence of the crossover scale on |ti |L is a signature of the one-dimensional 
nature of the large scale fluctuations in the fllm (recall In^ oc in a classical Ising chain) 
P|: long range order in a strip of flnite width is lost due to activation of two-dimensional 
domain walls running across the fllm; the Boltzman factor corresponding to creation of a 
pair of these walls is exp[—2Il{T)L/khT], where the free energy per unit length of the wall, 
S, happens to be given by in this formalism |TDI. 



I conjecture now that this mechanism of creation and annihilation of pairs of kinks at 
maxima and minima of g{z) deflnes the asymptotic RG flow for a generic temperature fleld 
t{z): at a given small q a solution can be constructed by matching pieces of (6) of different 
signs with the kink-solutions (7). A pair of kinks is stable while these kinks are separated by a 
potential difference Ag > | lng|; the pair is annihilated at smaller g's. This scenario correctly 
describes intermediate-scale behavior in the exactly solvable periodic LPI superlattice ||10|| . 
Because the fluctuations of potential g in the periodic model are strictly bounded, the critical 
flow is ultimately attracted by the uniform critical flxed point (4). On the other hand, a 
random, short range correlated t{z) gives rise to unbounded fluctuations of the random walk 



type: Ag'^(z) = D\z\ + o{z), with overbar standing for disorder average. Consequently, at 
a given q the typical size of a pair is going to be Az = Ingp, implying anisotropic 
scaling oc exp{D^/^^y^) at the critical random flxed point; the subscripts indicate 
directions parallel and perpendicular to the layers. Consider now a random multilayer grown 
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by susequent deposition of two different magnetic materials: t{z) takes two values, ti and 
^2, ti < t2 with probabilities p and 1 — p, respectively. The average t = pti + (1 — p)t2 goes 
through zero at a certain temperature T^, elsewhere a linear bias iz has to be added on top 
of the random walk exhibited by by g{z) at Tc. Crossover from the critical random to one 
of the noncritical fixed points happens when potential difference Ag{C,±) = iC,± coming from 
the linear part of g exceeds the typical random fluctuation on that scale, {D^±Y^'^ ~ In^y, 
so that on larger scales g is mostly monotonic. This consideration gives 

a = Di-^ C\\=exp{Dt-'), (11) 

in agreement with the identification of the average correlation lengths in [1(a)]; these lengths 
are different from the typical correlation lengths, as disorder-average of the logarithm of 
the correlation function results in ^_|_ oc i^^ Although, for instance, at t > the 



overall positive slope of g dominates at scales larger than ^, sufficiently long sequences of 
n > |ln(g)|/|ti| layers of pure low-temperature component occur with small but nonzero 
density p". The associated rare fluctuations in the slope of the potential keep a dilute 
system of large scale kink pairs from annihilation at any q. First appearance of these large 
scale pairs at the critical temperature of the first component, where ti becomes negative, 
leads to a Griffiths singularity in the free energy, which is easily estimated in agreement with 
0]; naturally, the same argument yields a singularity at Tc2- One should note that while 
the singularity at depends only on the universal large-scale variation, Ag oc z^^^, the 
essential singularities are governed by the tails of the distribution; if the distribution of t{z) 
is Gaussian the singularities are pushed all the way to T = 0, oo. 

This picture is easily generalized to the cases of long-range correlated random and 
quasiperiodic LPI models. In the first case, the main difference comes with the anoma- 
lous scaling dimension Xg > ^ describing the fluctuations Ag{z) oc z^^ at criticality, 
i = 0. The correlation lengths then diverge as oc t^"^, oc exp [const t^'^"], with 
i/_L = 1/(1 — Xb), t'li = Xb/{1 — Xh). On the other hand, if a quasiperiodic one-dimensional 
binary lattice is constructed by projecting from a strip of finite width cut out of a two- 
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dimensional lattice |TT|, then g{z) can be obtained by a similar projection and turns out to 
be confined within a strip of finite width around the average g oc iz. Thus critical behavior 
in this "Fibonacci" multilayer is governed by the pure Ising fixed point (4). Note that the 
amplitudes of the long wave length harmonics tkexp{ikz) of t{z) in the Fibonacci sequence 
are bounded by = in order for a quasiperiodic sequence to extrapolate between pe- 

riodic and random behavior, like it was found in for a quantum quasiperiodic XY-chain, 
the amplitudes have to go to zero slower than | A; | so as to generate unbounded fluctuations 
in the potential g. 

In conclusion, an exact functional real space renormalization group flow has been set 
up and studied for general layered planar Ising models. The RG scheme allows to separate 
relevant and irrelevant features of a wide class of systems, leading to the classification of 
various types of critical behavior based on simple basic characteristics, such as the asymptotic 
variation of the thermal potential g = 2 J tdz. The surface and random geometries have been 
analyzed within a single framework leading to a new physical picture of the criticality in 
a general disordered LPI model: At any finite momentum transfer q along the layers the 
system behaves as an alternating sequence of domains of the high- and the low-temperature 
phases. Having only one infinite dimension, the domains undergo an infinite coarsening as 
q decreases. The mechanism of coarsening is activation of linear defects across the finite 
widths of the domains, hence the exponential scaling (11) at the random criticality. The 
large fluctuations in the bond strengths underlying the domain structure in this formalism, 
must be essentially similar to those playing central role in the RG approach of [1(a)]. 

Discussions with Daniel Fisher, Mark Oxborrow, and Heiko Rieger have been very help- 
ful. 
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